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We study the ground state properties of the bond alternating S = 1/2 quantum 
spin chain whose Hamiltonian is 

H = + S y 2j Sl J+1 + \S* 2j S* w ) + 0Y1 S:y-i • S 2j • (0.1) 
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When (3 = 0, the ground state is a collection of local singlets with a finite excitation 
gap. In the limit of strong ferromagnetic coupling (3 — > — oo, this is equivalent to 
the S = 1 XXZ Hamiltonian. It has several ground state phases in the \-(3 plane 
including the gapful Haldane phase. They are characterized by a full breakdown, 
partial breakdowns and a non-breakdown of the hidden discrete Z<i x Z2 symmetry. 
The ground state phase diagram is obtained by series expansions. 
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I. INTRODUCTION 



We study the ground state properties of an S = 1/2 alternating spin chain with the 
Hamiltonian 

H = XX^j^+i + ^fjSfj+i + ^IjSajj+i) + P^2S 2 j-i ■ S 2 j, (1.1) 
j j 

where Sj = (Sj, Sj, Sj) are 5 = 1/2 spin operators. This model has two exchange couplings 
1 and (3 alternately. When (3 = 1 and A = 1, it is the spherically symmetric S = 1/2 
Heisenberg model which is solvable by the Bethe ansatz [jlj. When (3 = 0, the ground state 
is composed of local singlets that are formed by spins on sites 2j and 2j + 1. In the limit 
of strong ferromagnetic coupling (3 — > — oo, spins on sites 2j — 1 and 2j form a local triplet. 
The first order degenerate perturbation in 1/(3 from this limit gives the S = 1 spin chain 
with the Hamiltonian 

H 3 ' 1 = S*Si +1 + SfSf +1 + \SfS* +1 ), (1.2) 



where Sf, Sf and 5? are 5 = 1 spin operators at site i. When A = 1, the Hamiltonian ( |1.2| ) 
is the Heisenberg antiferromagnetic chain with S = 1. 

In 1983 Haldane conjectured that there are qualitative differences between integer and 
half-integer spin quantum antiferromagnetic chains based on the large-S" expansion [Q. 
While there is no excitation gap when S is a half- integer he argued that the Heisenberg 
Hamiltonian has a unique disordered ground state with a finite excitation gap when the 
spin S is an integer. These claims have been confirmed by experimental [§J , numerical 
and analytical studies at least for the S = 1 case. The elementary excitations and the 
temperature dependence of the Haldane phase are investigated using the defect approach 
|T0| . Although the ground state is disordered in the sense that the usual spin-spin correla- 
tion function decays exponentially, it has turned out that the ground state has a nontrivial 



hidden order. Den Nijs and Rommelse ['I I ] argued that the Haldane phase is characterized 



by the hidden antiferromagnetic order based on the analogy of the preroughning transition. 
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It is measured by the string order parameter. Kennedy and Tasaki [|12j introduced a nonlo- 
cal unitary transformation which reveals the relation between the hidden antiferromagnetic 
order and the hidden discrete Z2 x Z2 symmetry breaking. The hidden Z2 x Z2 symmetry 
is completely broken in the Haldane phase. This is confirmed numerically |13[. Extensions 
to higher integer S cases were discussed by several authors (H] . 

The S = 1/2 spin chains with bond alternation have been investigated by several groups 
1 15 1 - Recently Hida studied an S = 1/2 Heisenberg chain with alternating ferromagnetic 



and antiferromagnetic couplings. Its Hamiltonian is ( [LI ) with A = 1. Hida considered 
a string order parameter for S = 1/2 which leads to the den Nijs-Rommelse string order 
parameter in the (3 — > — 00 limit. He concluded numerically that the gap and the string 
order parameter for (3 = remain non- vanishing in the (3 — > — 00 limit. It suggests that 
the S = 1/2 spin chains with disordered ground state are characterized by the string order 
parameters, and that for A = 1 there is no phase transition between f3 = (S = 1/2 local 
singlet-triplet gap) and (3 = —00 (S — 1 Haldane gap), thus they belongs to the same phase 



There are studies on the 5 = 1/2 spin chains with bond alternation using other types of 
nonlocal unitary transformation JT9|,pD| . In particular, Kohmoto and Tasaki [T9j studied the 



Hamiltonian (|1.1|) by the nonlocal unitary transformation used in the work of the Ashkin- 
Teller model by Kohmoto, den Nijs and Kadanoff |21]]. They found that this transformation 
plays a role of the Kennedy- Tasaki unitary transformation for the S = 1 chain, although 
they are distinct. The nonlocal unitary transformation maps the S = 1/2 alternating chain 
to an S = 1/2 quantum spin system which is similar to the highly anisotropic version (one- 
dimensional quantum system) of the two-dimensional Ashkin- Teller model. It also maps the 
string order parameters of the original chain to the ferromagnetic "local" order parameters 
in the transformed system. They argued the correspondence between the breaking of the 
Z2 x Z2 symmetry in the transformed system and the hidden antiferromagnetic order in 
the original system. In the transformed system, they showed that the region including the 
decoupled model ((3 = 0) and the S = 1 chain are characterized by a full breakdown of the 



Zi x Z2 symmetry. The Neel ordered phase (in terms of S = 1) is characterized by a partial 
breakdown. They also proposed a ground state phase diagram. 

The unitary transformation maps the quantum spin chain with bond alternation into the 
Ashkin- Teller type quantum spin system on a pair of chains. The spin chains are similar to 
the Ashkin- Teller quantum chain. The latter model is obtained by a highly anisotropic limit 
of the two-dimensional Ashkin- Teller model The two-dimensional Ashkin- Teller 

model |23| consists of two Ising models on a square lattice coupled by a four-spin inter- 
action. In the time-continuum Hamiltonian formalism, a two-dimensional classical system 
is reduced to a one-dimensional quantum system by taking the extreme lattice anisotropic 
limit. The transfer- matrix method is used to convert a statistical mechanics problem at a 
finite temperature in two dimensions into a ground state problem for one-dimensional quan- 
tum Hamiltonian p4| . In this formalism, the ground state energy of a quantum system is 
the free energy of the corresponding two-dimensional classical system. The excitation gap 
in the quantum system corresponds to the inverse of the correlation length of the classical 
system. If there is no excitation gap in the quantum system, the classical system is critical. 
Namely phase boundaries of the ground states phase diagram of a one-dimensional quantum 
system correspond to critical points of the two-dimensional classical spin system. Thus we 
shall use the concepts of the critical phenomena (universality, critical indices etc.) in order 
to discuss the ground state phase diagram of the one- dimensional quantum system even in 
the parameter region where a corresponding two-dimensional classical system does not exist. 

The purpose of the present paper is to determine the phase diagram of the Hamiltonian 
( |1.1|) quantitatively by a series expansion technique. Some limiting cases are considered 
analytically. The critical lines and critical indices are evaluated from the Pade method ||25| . 



The obtained phase diagram shows a rich structure. It contains the gapful Haldane phase, 
the Neel phases, the XF-like gapless phase. It also has a critical line with continuously 
varying critical indices. 
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II. THE DISORDER OPERATORS AND THE Z 2 x Z 2 SYMMETRY BREAKINGS 



We shall perform series expansions in terms of (3. The unperturbed system ( ( |1 . 1| ) with 
(3 = 0) has the ground state which is a collection of uncorrelated local singlets. It is 
disordered in the sense that the expectation values of local order parameters vanish and 
series expansions cannot be applied to these quantities. Thus we consider the disorder 
operators which are nonzero in disordered phases and zero in an ordered phase We 
choose the following disorder operators 

D XV {j) = Y\{ —^(SZkSZk+l + ^2k^2k+l) } = II { ~(^2k^2k+l + ^2k^2k+l) }> (2-1) 
k>j k>j 

and 

D z (j) = U(-^ 2k S^ +1 ), (2.2) 

k>j 

to characterize the ground states of the Hamiltonian ( |1 . 1|) , where S a are the 5 = 1/2 spin 
operators. 

We describe the derivation of above disorder operators and their relations to the hidden 
Z2 x Z2 symmetry. The nonlocal unitary transformation reveals the Z2 x Z2 symmetry of 
the Hamiltonian ( | . 1|) . This transformation is exactly the same as that used by Kohmoto, 
den Nijs and Kadanoff jnj], which maps the staggered XX Z model into the Ashkin- Teller 
quantum chain. Applying this transformation to the Hamiltonian ([Lip, we get 



H = ~ + ^ + A*?7f ) - P £(^f + i + t*t*+i + ^< + i<rf +1 ), (2.3) 

i i 

where af and r" (a = x, y and z) are different kinds of Pauli matrices at site i. For details, 
see Appendix B of Ref. 21. 

This Hamiltonian has symmetries of the dihedral group of order 4 which includes several 
Z2 symmetries. There are many varieties of order parameters to measure the spontaneous 
symmetry breaking of these symmetries. We choose the following order operators: 

0± = -s-p-, (2.4) 
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2 = olrl. (2.5) 



The disorder operators which are dual to these order operators are 



D ± (i) = U^ 11 , (2-6) 

k>i A 

Ai(i) = II^ ( 2 - 7 ) 

k>i 

respectively. These disorder operators have nonzero expectation values in disordered phases. 
Applying the inverse of the above transformation to the disorder operators Q2.6|) and 
we get the disorder operators (|2.1|) and ( |2.2|) of the Hamiltonian ( |1 . 1|) . 



A nonzero expectation value of one of the disorder operators corresponds to a spontaneous 
breaking of a Z 2 symmetry. The ground state phases are characterized by a full breaking, 
partial breakings and nonbreakings of the Z 2 x Z 2 symmetry. Thus the ground state phase 
diagram can be obtained from analysis of the disorder operators ( |2.1| ) and (|2.2| ). 

The disorder operators ( j2TT|) and ( |2.2|) are related to the string order parameters of den 
Nijs and Rommelse, and Hida. It is expected that they have the same critical properties. 

III. PHASE DIAGRAM 

Before going to the detailed analysis, we summarize our results by showing the obtained 
phase diagram in Fig. |I[ This phase diagram have seven phases and eight critical lines. 

(1) The gapful Haldane phase (A). It is characterized by ( D a ) 7^ 0, (a = xy or z). 

(2) The S = 1 Neel phase (B) which is characterized by ( D z ) 7^ and ( D xy ) = 0. This 



region continuously connects to the Neel phase of the S = 1 Hamiltonian ( |1.2|) . 

(3) The AF-like gapless phase (C) where ( D a ) (a = xy and z) vanishes. The correlation 
function decays algebraically in the entire region and the critical indices vary continuously. 

(4) The ferromagnetic phase (D). 

(5) The disordered phase (E) in which the ground state is a disordered dimer state with 
( D a ) =0 (a = xy and z). 
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(6) The S = 1/2 Neel ordered phase (F). This phase is characterized by ( D z ) ^ and 
(D*v ) = 0. 

(7) The Neel ordered phase (G) where spins on sites 2j and 2j + 1 prefer to align parallel 
and those on sites 2j — 1 and 2j prefer to align antiparallel. 

These phases are separated by the following critical lines. 
(1) Line 1 is expected to be in the Gaussian model universality class |27|. Critical indices 



vary continuously on this line for — 1 < A < 1. This critical line is expected to bifurcate at 
A = l§. 

(2) Line 2 and Line 3 are expected to be in the Ising model universality class. Line 2 ap- 
proaches the line (3 = 2 for large A. Line 3 approaches the line f3 = |A for large A. 

(3) Line 4 is the phase boundary between the Haldane phase and the Neel phase (in terms 
of the S = 1 model) and belongs to the Ising model universality class. 

(4) Line 5 is the phase boundary between the XF-like gapless phase and the Haldane phase, 
and is expected to be in the universality class of the Kosterlitz-Thouless transition ||29|| . 

(5) Line 6 (A = — 1, (3 < 0) is the phase boundary between the XF-like phase and the 
ferromagnetic phase, and is expected to be in the universality class of the KDP transition. 



(6) The Hamiltonian ( |1 . 1|) is solvable on line (3 = 0. The ground state is degenerate on Line 



7. Thus, Line 7 is a critical line 21 



(7) Line 8 belongs to the Ising model universality class and approaches the line (3 = — |A, 
for large |A|. 



A. Limiting Cases 

In this subsection we discuss the phase diagram in several limits. Some of the critical 
lines are obtained by the duality. 
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1. limit A > \(3\ ~ 1 



Since the term with A dominates, we can restrict ourselves to the space spanned by 

I ) = <8> \ a i >2j,2j+l , (3.1) 

3 

where the state \aj) 2 j,2j+i is defined by 

I OCj =T )2j,2j+l = | T)2j | V)2j+1, | OCj =1 )2j,2j+l = | i>2j | thj+l- (3.2) 

Here, | f)j and | are eigenstates of 5| with eigenvalues +1/2 and —1/2 respectively. We 
treat the other parts of the Hamiltonian by the degenerate perturbation theory within the 
space spanned by fl3.1| ). The first order perturbation gives the effective Hamiltonian 

# = ^£<-~/3£^K +1 , (3-3) 
where af are the Pauli matrices which operate on \oti This is the highly anisotropic 



version of the two-dimensional Ising model ||24| . This Hamiltonian have the property 

2 \!3 



~^h(- o , ^)=H((3,a a ) , (3.4) 



under the dual transformation 

< = IT af = aM +l . (3.5) 

k=l 

For /3 < 0, this duality of the Ising model tells us that (3 = —2 is the critical line. Thus the 
critical line between the Haldane phase and Neel phase approaches (3 = —2 as A becomes 
large. Above this line (phase A in Fig. |l|), the Z 2 X Z 2 symmetry is fully broken. Below 
this line (phase B in Fig. |l[), the ground state breaks half of the Z 2 x Z 2 symmetry. Critical 
properties corresponding to a partial breakdown of Z 2 x Z 2 symmetry is considered as the 
Ising model universality class. Thus Line 4 belongs to the Ising model universality class. 

For (3 > 0, we apply a transformation S 2 j — > —S 2 j, Sfj "~ > — &2j an d S 2 j — > S 2 j only 
for spins on the 2j sites. Then we get the same duality relation as ( |3.4| ) except for the sign 
of (3. It implies that (3 = 2 is a self-dual line for the effective Hamiltonian. Thus Line 2 
which approaches (3 = 2 as A becomes large belongs to the Ising model universality class. 
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2. limit (3, A > 1 



The Hamiltonian is approximated by 

H = A E S^ +1 S*-i • S 2J . (3.6) 

When A is sufficiently large, the A-coupling dominates in the Hamiltonian ( |3.6| ) and spins on 
sites 2j and 2j + l must form a state | f ) 2j - 1 Ihj+i or \ [) 2 j\ T)2j+i- In this bases, spins on sites 
2j — 1 and 2j favors the states, | T)2j-i| D23 or I l)2j-i| T ) 2j , due to the antiferromagnetic 
coupling (3. Therefore, the ground state is given either by 1 \) 2 j | |)2j+i) or <8>j(| l)2j I T 
) 2j '+i). It is a doubly degenerate Neel ordered state. The energy per site is |A — On the 
other hand, spins on sites 2j — 1 and 2j favor to form singlet for f3 ^> A. The ground state 
is given by an array of them. It is disordered and the energy is — 1/3. 

To investigate the phase boundary between these phases, we map the Hamiltonian (|3.6|) 
into the model without bond alternation. Applying the nonlocal unitary transformation (see 
Sec. 0) to the Hamiltonian fl3.6|), we get (see (|2.3|) ) 

H = -\J2°M-f3 EK< + i + if rf +1 + « +1 r/rf +1 ). (3.7) 

i i 

We map the space spanned by the eigenstates of of and rf to that spanned by the eigenstates 
of ofif and if. We define the operator rf by rf = of-rf . Using this notation, the Hamiltonian 
( ft.Tp becomes 



H = ~^J2 T i ~ $ E(?Cl + ^1 + 7f 7f 7f +1 7f +1 ) (3.8) 

i i 

= -A E 7f - ^E{(! + + Ti-rUY (3-9) 

i i 

In this Hamiltonian, the coupling of rf is negative, since 1 + t't? +1 is always positive. Thus 
i spins appear as the one- dimensional Ising model with the ferromagnetic coupling. In the 
ground state, f spins are completely ordered ferromagnetically. We can replace the operator 
ifii+i by its expectation value (rfr/ +1 ) = 1. Thus we have 

iy = -AE^-2/3E^/+i> ( 3 - 10 ) 



apart from a trivial additive constant. The effective Hamiltonian ( |3.10| ) is the Ising model 
in a transverse magnetic field. Therefore, the boundary of the phases is A = |/3 in this limit 
(Line 3 in Fig [I]). Note that the phase transition for r spins across critical line 3. Since 
the Hamiltonian ( |3.7| ) is symmetric with respect to r spins and a spins, a spins have the 
same behavior as r spins. Therefore, the order parameter ( |2.4j ) has nonzero expectation 
value in Phase E. The disorder operator whose series has singularity on Line 3 is D z . Below 
this line (phase F in Fig |l|), the system is ordered in terms of S — 1/2 Neel order and is 
characterized by ( D z ) ^ and ( D xy ) = 0. Thus the ground state breaks half of the 
Z 2 x Z 2 symmetry. Above this line (phase E in Fig |l|) , the system is disordered and is 
characterized by ( D a ) = (a = xy and z). The Z 2 x Z 2 symmetry is fully restored. Line 
3 belongs to the Ising model universality class, since the Z 2 symmetry is broken across this 
line. 



3. limit P, |A| (A < 0) » 1 
Similar to the previous subsection, the Hamiltonian of this regime is 

h = -\MY: s* 2j si j+1 + ^Vi • S 2J • (3.11) 

3 j 

In the case of large |A|, spins on sites 2j and 2j + 1 prefer to form a state | t) 2 j| thj+i 
or | i) 2 j\ Dij+i- Spins on sites 2j - 1 and 2j favor the state | 1) 2j \ l) 2 j+i or | l) 2j \ f 
)2j+i, since the coupling j3 is antiferromagnetic. Therefore, the ground state is given by 
0j ( | T>4j I T>4j+i I i)ij+2 I i)4j+s ) and ( | [) Aj \ \ ]) ij+2 \ t>4y+3 )• The energy 

of this state is — ||A| — \(5 per site. For (3 3> |A|, spins on sites 2j — 1 and 2j favor a singlet 
state. The ground state is given by a collection of local singlet with the energy — 1/3. In this 
case, the phase boundary is the line A = — 1/3 (Line 8 in Fig. |I|). Below this line (Phase G 
in Fig. |I|), the system breaks half of the Z 2 x Z 2 symmetry. Thus the transition across this 
line belongs to the Ising model universality class. 
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4- line (3 = 



We have the simplest situation. The Hamiltonian is written 

H = /XS%jS% i+ -t + 5 , | J -5 , 2 J+1 + ^S^j^Zj+i) ■ (3.12) 
j 

It is a collection of independent two-spin systems. For A > — 1, the ground state is given 
by (I T>2j| l>2i+i - I V)2j\ T>2j+i)} with an energy per site. For A < -1, the 

ground state of an independent two-spin system is given by | ]) 2 j\ t)2j+i or | [) 2 j\ D^j+i- 
Therefore, the ground state of the total system is 2 M//2 fold degenerate, where M is the 
number of sites. Therefore, this is a critical line [EI]]. The energy is | per site. 



5. the ferromagnetic phase (3 < 0, A < — 1 

In this region, the Hamiltonian is written 

H = XX ^2? £27+1 + S 2j S 2j+1 - \M^2j^ij+i) - ^2 ^2j-i • S 2j f • (3.13) 

j 3 

The ground state is given by 3 - 1 1) 2 j\ t)2j+i or 0^ | l) 2 j\ D2J+1 with the energy - |A| g l/31 . 
Line 6 is expected to be a critical line in the universality class of the KDP transition 



3tJ, since the system is in the perfectly ferromagnetic ordered state in the A < — 1 side and 



is in the XF-like gapless phase in the A > — 1 side. 



6. limit (3 — > — 00 

In the (3 — > —00 limit, the S = 1/2 Hamiltonian (|1 . 1| ) is equivalent to the S = 1 XX Z 
Hamiltonian ( |1 . 2| ) . It has the Neel, the Haldane, the XF-like and the ferromagnetic phases 
depending on the parameter A. We denote Ai for the boundary between the Neel and the 
Haldane phases, A2 for that between the Haldane and the XF-like phases, and A3 for the 
XY-\ike and the ferromagnetic phases. 

The phase boundary between the Neel and the Haldane phases is estimated to be A 1 ~ 1.2 
numerically The critical line line 4 is expected to approach A = X\ [[0|]. The boundary 
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between the Haldane and the XY-\ike phase is estimated A2 ~ numerically ||||. The 
XY-like phase boundary which starts from the point (A, 0) = (—1,0) approaches A = A 2 . 

For A = — 1, the 5=1 XX Z Hamiltonian ( |1.2j ) can be mapped to the ferromagnetic 
Heisenberg model by the rotation in the spin space for the 2j sites, 5fj- — ► —5^, 5|j — > — 5fj 
and 5|j — »■ 5f-. The ground state have degeneracy due to the rotational symmetry at 
A = — 1. On the other hand, for A < —1, the ground state is 2 fold degenerate. Therefore, 
the phase boundary between the XY-\ike and the ferromagnetic phases is A = —1 in this 
limit. 



IV. SERIES EXPANSION 

We will make a series expansion in terms of (3 and obtain series of the "specific heat", 
the "magnetization" and the "susceptibility" to estimate critical points and critical indices. 



A. Outline of the Method 

It is necessary to have long series to extract reliable estimates of the critical points and 
critical indices. We use the linked cluster expansion method proposed by Kadanoff and 



Kohmoto It is suited to calculate high order terms effectively. 



Hamiltonian fll.lQ is written 

H = T + (3V, (4.1) 



where 



T — / ] (51,51,4.1 + ^2j^2j+i + ^2j^2j+i) 5 (4-2) 

3 



and 



K = X)S 2 i-i-S 2i . (4.3) 



j 



We shall perform series expansions with respect to (3. The unperturbed states are 
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W ) = <S> Wjhj, 23+1, (4.4) 



where \oij)2j,2j+i is one of four eigenstates of the operator S 2j S 2j+1 + S 2j S 2j+1 + ^S 2j S 2j+1 : 

|'S')2j,2j+l = 75 (| T ) 2 j | Ihj+l - I l)2j| T)2j+l), |7+l)2j,2j+l = | T > 2jr | t^j+l, |^o)2j,2j+l = 

75 (I T>2j| i>2i+i + I !>2j| T>2j+i) and \T_ 1 ) 2j}2j+1 = \ l) 2j \ [) 2j+1 . For A > -1, the unper- 
turbed ground state is 

\G) = <g)\S) 2] , 2H1 . (4.5) 

3 

The ground state energy is expanded in a power series as 

E(\,P) = Y / E (n) (\)p n . (4.6) 

?1=0 

The "specific heat" is obtained by 



C{\(3) = ^E{\(3). (4.7) 



We use the disorder operators defined per site 



^ = ^£^), ( 4 - 8 ) 

3 

where a = xy and z, and M is the number of sites. To calculate series for the "magnetiza- 
tion" of the operators, we simply add a magnetic field 

H = T + (3V + hV a . (4.9) 

The ground state energy in the presence of h is calculated in a power series of (5 and h as 

E(X;P,h) = E E E {n ' m \\)f3 n h m . (4.10) 

n=0 m=0 

The expectation value of the "magnetization" is given by 

m = , (4.11) 

and the "susceptibility" is 
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d{ V a ) d 2 E(X; {3, h) 



dh 



dh 2 



\h=0 



(4.12) 



We use the Pade method to estimate critical points (3 C and critical indices. The definition 
of critical indices are as follows: 



8F_ 

dp 



-E(X, (3) ~ constant 



Pc~P 



Pc 



(P^Pc), 



(4.13) 



dE{\; p, h) 
dh 



\h=0 



constant 



Pc-P 



Pc 



(P^Pc), 



(4.14) 



d 2 E(\;P , h) ] 



\h=0 



constant 



Pc-P 



Pc 



~7j 



(/?-&), 



(4.15) 



where j is xy and z. 



B. Result from Series Analysis 

The quantities calculated by series expansions were (1) the "specific heat" (13 terms), 
(2) the "magnetization" (15 terms), (3) the "susceptibility" (13 terms). Critical points 
and critical indices are evaluated by the Dlog Pade method ||32|| . The estimates for the 
quantities are obtained by averaging the three or four highest-order diagonal elements and 
near-diagonal elements of the Pade tables. Error bars are then set to include these three or 
four values. 



1. region (3 > 

For the region — 1 < A < 1, the Haldane phase is characterized by (D a ) ^ (a = xy 
and z). The disordered phase (E) is characterized by (V a ) = (a = xy and z). Thus the 
phase boundary (Line 1) between these phases is determined from series for T> a (a = xy 
and z). Critical points obtained from series for (T> xy ) and (T> z ) show well convergence. We 
estimate critical points from series for (V xy ) (Fig. [I]). In Fig. and|3|, critical indices P/s (j 



14 



= xy and z) are shown. Results for 7j's (j = xy and z) and a are shown in Fig. f|, |5] and |6| 
respectively. In this region, critical indices /3-,'s (j = xy and z) vary continuously, and seem 
to be divergent near A = —1. Critical index j3 z agrees well with the result by Hida obtained 
numerically |[33|| . Critical indices 7/s (j = xy and z) and a are also varying continuously. 
Critical index 7^ also seems to be divergent near A = — 1. However, we cannot determine 
critical indices accurately close to A = — 1 due to poor convergence. In the case of the 
Ashkin- Teller quantum chain or the staggered XXZ model, there is the "critical fan" in 
a finite region of the parameter space — 1 < A < — -j= ^T\. The "critical fan" is a region 
where a line of continuously varying criticality "fans out" and becomes an area of critical 
behavior. In that case, critical indices show divergence near A = — A= on the self-dual line of 
the Ashkin- Teller quantum chain. In the present model, there is no sign of a "critical fan" 
within this analysis, since critical points show well convergence and critical indices does not 
show divergence. However, we cannot discuss the existence of the "critical fan" near A = — 1 
for poor convergence of critical indices. 

For (3 = 1, A = 1, the model is solvable by the Bethe ansatz We can map it to the 
six- vertex model at this point ||34|| . The extended scaling relations predict critical indices 
a, /3j's and j/s (j = xy and z), with the help of the standard scaling relations. These 
indices have been obtained from mappings between the six-vertex model and the Gaussian 



model p7|j35fl . The equivalent mapping is that between the XXZ model and the Tomonaga- 
Luttinger model [j36|| . The estimated values of the critical point from the series analysis are, 
P c = 1.003 ± 0.002 (from series for the "specific heat"), (3 C = 1.00168 ± 0.00004 (from series 
for (V xy )), (3 C = 0.99782 ± 0.00009 (from series for (V xy V xy ) ). The critical indices a, f3/s 
and Tj's (j = xy and z) are 

a = 0.55 ±0.01 (^= 0.666--), (4.16) 

Pj = 0.08202 ± 0.00008 (— = 0.0833 ■ •■ ) , (4.17) 

7 

7 3 - = 1.23 ±0.02 (- = 1.166--) ■ (4.18) 
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The values in parentheses are the expected values from the mappings. 

Next we consider the region A > 1. The Haldane phase is characterized by (V a ) ^ 
(a = xy and z). The 5 = 1/2 Neel phase (F) is characterized by (T> xy ) = and (T> z ) ^ 0. 
Thus the phase boundary (Line 2) between these phases is determined from series for T> xy . 
In Phase E, both of (V xy ) and (T> z ) vanish. Therefore, the phase boundary (Line 3) is 
determined from series for D z . These results are shown in Fig. |l|. Line 2 and Line 3 belong 



to the Ising model universality class (see Sec. |T|). Thus critical indices a, (3 xy and ^ xy must 
be 0, | and |, respectively on Line 2. Critical indices a, (3 Z and j z must be 0, | and |, 
respectively on Line 3. The series analysis shows that Line 2 approaches the line j3 c = 2, 
which is consistent with the result in Sec. |HI[ For A > 1, Critical index f} xy agrees well 
with that of the Ising model |. Critical index ^ xy also agrees well with that of the Ising 
model j. Near A = 1 they do not agree with the Ising values. We regard this as a numerical 
effect. Critical index a does not agree well with the Ising value a = as A becomes large. 
In the series analysis for (T> z ), there is a pole near Line 2. We choose the second pole to 
determine the critical points and critical indices for Line 3 as shown in Fig. [I]. The operator 
T> z has nonzero expectation value in Phase A and in Phase F. We regard the first pole as 
unphysical consequences, since series for T> z has no singularity on Line 2. The series analysis 
shows that Line 3 approaches the line /3 = |A. This is consistent with the result in Sec. 



nj. However, their convergence is poor. Critical indices a, (3 Z and 7 2 are not included here 
for poor convergence. Critical indices f3 z estimated from series analysis are f3 z = 0.03 ~ 0.1 
for A = 2 ~ 5. They do not agree well with the Ising values. Critical index ^ z does not 
converge. We think that these are due to the first pole near Line 2 and large value of (3 C . 



2. region (3 < 

Let us consider the region A > 0. The 5 = 1 Neel phase (B) is characterized by (T> xy ) = 
and (T> z ) ^ 0. Thus the phase boundary (Line 4) between the Haldane and the 5 = 1 Neel 
phases is determined from series for D xy . The best estimates are obtained from series for 
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(T> xy ) as shown in Fig. |l|. Line 4 belongs to the Ising model universality class (see Sec. 
Therefore, critical indices are expected to take the Ising values throughout Line 4. In the 
limit (3 — > — oo, Line 4 is expected to approach the Neel-Haldane transition point Ai of the 
S = 1 XXZ model. The critical points and critical indices obtained from series for (V xy ) 
are well convergent for A ~ 2.5, and stable for 2 ~ A ~ 2.5. However, they are unstable for 
A ~ 2. No pole was found in every highest-order element of Pade tables for < A ~ 1.2 
except for them with small residue (< 10~ 5 ). Note that if we investigate the behavior of 
Pade tables carefully, they show a highly unstableness near A ~ 1.8. Thus, for the region 
1.2 ~ A ~ 1.8, we regard the critical points as poor convergence of the Pade approximation. 
For large A, the critical line approaches the line f3 c = —2, which is consistent with the result 
in Sec. |T|. Critical index (3 xy is well convergent to that of the Ising model | (see Fig. 0). 
The result for ^ xy is shown in Fig. [8[ It shows a good agreement with the Ising value %. 
Convergence of critical indices become poor around A = 2. This is due to large value of /3 C . 

Next we focus on the region — 1 < A < 0. In the XF-like phase, both of (T> xy ) and 
(T> z ) vanish. Thus the phase boundary (Line 5) between the Haldane and the XF-like 
phases is determined from series for T> a (a = xy and z). We choose the operator D z to 
determine Line 5, since critical points obtained from it converges better. The critical points 
and indices are stable near the point (A, (3) = (—1, 0). The convergence of them becomes 
poor when A approaches 0, and it is difficult to determine the critical points and critical 
indices. Continuously varying critical exponents are observed for (3 xy and (3 Z along this line, 
which are shown in Fig. [7| and Fig. |9| respectively. Critical indices a and 7-,'s (j = xy and 
z) are not included here for poor convergence. These values estimated from series analysis 
are a = 2.2 ± 1.1, ^ m = 3.4 ± 2.6 and 7 2 = 3.1 ± 2.6 at A = —0.8. Line 5 is expected to be 
in the universality class of the Kosterlitz-Thouless transition. However, we cannot evaluate 
critical index rj for poor convergence of critical indices jj's (j = xy and z). 
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V. SUMMARY AND DISCUSSION 



We studied the ground state properties of the 5=1/2 quantum spin chain with bond 
alternation. This model is equivalent to the 5 = 1 antiferromagnetic XXZ model in the 
strong ferromagnetic coupling limit. We focused on the disorder operators. These operators 
are also obtained from the disorder operators of the Ashkin- Teller quantum chain by using 
the nonlocal unitary transformation. The "magnetization" of these nonlocal operators mea- 
sure similar quantity that is observed by the string order parameters which is extended by 
Hida for the 5=1/2 quantum spin chains. 

The ground state phase diagram is obtained by series expansions. Critical points are 
consistent with the phase diagram which conjectured by Kohmoto and Tasaki and 



determined numerically by Hida ||18|| . In the narrow region A > — 1, (3 ~ 0, it can be proved 
rigorously that the ground state in the infinite volume limit is unique, all the truncated 
correlation functions decay exponentially, and there is a finite excitation gap |12| , |l9| . The 
present study strongly suggest that this region covers the wide area surrounding by Lines 1, 
2, 4 and 5 in Fig. [|. The results of critical indices show that phase transition between the 
Haldane phase and the Neel phase (Line 4 in Fig. [I]) belongs to the Ising model universality 
class. The Ising-like critical properties are originated from the partial breakdown of the 
Z 2 x Z 2 symmetry in this case. This supports that the hidden Z 2 x Z 2 symmetry breaking 
is to be a criterion to distinguish the Haldane gap system. Although the nonlocal unitary 
transformation discussed here is not exactly the same that used by Kohmoto and Tasaki, 
the disorder operators have same critical properties with theirs. This result also supports 
the proposal by Hida that the Haldane phase which is characterized by the string order 
parameter belongs to the same phase as that of the decoupled 5 = 1/2 chain. 
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FIGURES 

FIG. 1. Phase diagram of the alternating quantum spin chain with the Hamiltonian fll.ip , 
Estimates of critical points by a series analysis are shown with error bars. For those without an 
error bar, the error is smaller than the size of the plotted point. 

FIG. 2. Critical index (3 xy as a function of A on Line 1 and 2. 

FIG. 3. Critical index f3 z as a function of A on Line 1. 
FIG. 4. Critical index j xy as a function of A on Line 1 and 2. 

FIG. 5. Critical index ^ z as a function of A on Line 1. 
FIG. 6. Critical index a as a function of A on Line 1 and 2. 
FIG. 7. Critical index (3 xy as a function of A on Line 4 and 5. 

FIG. 8. Critical index ^ xy as a function of A on Line 4. 

FIG. 9. Critical index (3 Z as a function of A on Line 5. 



23 



